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We propose a framework for calculating scattering and bound state properties in anisotropic two-dimensional 
potentials. Using our method, we derive systematic approximations of partial wave phase shifts and binding 
energies. Moreover, the method is suitable for efficient numerical computations. We calculate the s-wave phase 
shift and binding energy of polar molecules in two layers polarized by an external field along an arbitrary 
direction. We find that scattering depends strongly on their polarization direction and that absolute interlayer 
binding energies are larger than thermal energies at typical ultracold temperatures. 
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I. INTRODUCTION 

Anisotropic interactions are at the heart of many physical 
systems, such as the atomic nucleus, heteronuclear molecules, 
or ultracold dipolar atoms. The anisotropy has a marked in- 
fluences on their scattering and bound state properties. Ad- 
ditionally, confining a system to layers often leads to exotic 
effects and is thought to be a crucial ingredient for the elusive 
theory of high- T c superconductivity [1]. 

Heteronuclear molecules are a prime candidate for study- 
ing low-energy anisotropic scattering in a controlled manner. 
Notably, they possess a large permanent dipole moment. Fur- 
ther, several species have now been cooled to ultracold tem- 
peratures [2]. These conditions give rise to an anisotropic 
dipole-dipole interaction (DDI), which can be controlled by 
external fields [3]. For example, an electric field polarizes 
the molecules. However, if the DDI dominates the three- 
dimensional (3D) dynamics, the gas becomes unstable [4]. 
Confining the molecules to a two-dimensional (2D) layer sta- 
bilizes the gas for perpendicular polarization [5]. Perpendicu- 
larly polarized molecules in two parallel layers give rise to in- 
terlayer bound states and scattering across layers [6-8]. Vary- 
ing the polarization direction with an external field controls 
the anisotropy of their interaction. However, only a very re- 
cent work has studied the influence of this DDI anisotropy on 
bound states [9]. 

In this paper, we propose a general framework for calcu- 
lating scattering phase shifts and binding energies for a wide 
class of 2D potentials including anisotropic potentials and the 
special case of vanishing potential volume J d 2 rV (r) = 0. Our 
method and calculation extend the results in a central 2D po- 
tential by Klawunn et al. [6] to the anisotropic case. Bound 
states and scattering in other classes of 2D potentials have 
been studied in the past [10, 11]. On the one hand, our for- 
malism allows for efficient numerical computations in 2D as 
it only requires solving a first-order differential equation. On 
the other hand, we present systematic approximations for scat- 
tering phase shifts and binding energy in anisotropic 2D po- 
tentials. Specifically, we recover the Born approximation for 
the phase shifts and an approximation connecting the binding 
energy and the total phase shift at low energies. As an ex- 
ample, we calculate the interlayer binding energy and s-wave 
phase shift of polar molecules in a bilayer at arbitrary polar- 
ization. The interlayer binding, e.g., of two 6 Li 40 K molecules, 



is deeper than thermal energies in typical ultracold experi- 
ments over a range of polarization angles. 

We consider two particles interacting via a 2D poten- 
tial V(r). In the center-of-mass frame the scattered state 
depends on both the relative vector r = r(cos(p,sin(p) be- 
tween the two particles and the incoming wave vector k = 
fc(coscj,sinc^). Consequently, we expand the dimensionless 
wave function in partial waves for both vectors: ^(kjr) = 
2i £m n=-°° ^"'"l^^ e"'"Pe~ m ^. Inserting this expansion into the 
radial Schrodinger equation of the two particles results in 
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Here, k = yE2fjA 2 /h is the collision momentum at en- 
ergy E, /j is the reduced mass, and A is the unit of 
length. The matrix elements of the potential are V mm i(r) — 

Different physical solutions of Eq. (1), such as scattering 
or bound states, fulfill different boundary conditions. For our 
proposed formalism, we expand all physical states \|/ m „ in a set 
of regular solutions of Eq. (1), (]),„„, with well-defined bound- 
ary conditions. Specifically, we require that these regular so- 
lutions behave at the origin as lim,^o§ nm (k,r) / j n (kr) = 8 mn . 
Here, we have defined the scaled Bessel functions j m (x) = 
\Z%x/2J m (x) and y,„ (x) = yJnx/2Y m (x), where J m and Y m are 
Bessel functions of the first and second kind, respectively. For 
2D scattering, these scaled Bessel functions are the equiva- 
lent of the Riccati-Bessel functions familiar from 3D scatter- 
ing [12]. Our choice of the boundary condition fixes the free- 
dom in the expansion coefficients of the regular solution and 
its derivative [13]. 

We assume \im r ^o r 2 V m „(r) = and lim^oo rV mn (r) = 0. 
Because of the latter, at large distances solutions of Eq. (1) are 
proportional to the scaled Hankel functions h^(x) = j m (x) ± 
iy m (x). This motivates us to expand the regular solutions as 

4W, (k,r) = \ K (kr)ti„ (k,r) + h m (kr)f mn (k,r)]. (2) 

We insert this expansion into Eq. (1) and require that 
h+(kr)d r f+,(k,r)+h m (kr)d r f mn (k,r) = 0. This condition re- 
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duces the Schrodinger equation to the first order equation 

3 r /4(M = i^^VW^^r) (3) 



for the coefficients ff m . The formal solution of Eq. (3) is 
1 



fmn( k ^) = §m«± / dr' h+(kr')YV mm i (r')^ m > n (k,r'). 

(4) 

Here, we have fixed the integration constants so that f^ n (k, °°) 
exhibits the correct high-energy behavior [13]. 



II. BOUND STATES 

A weak 2D potential supports bound states if / d 2 rV(r) < 
[11]. Using the coefficients f^ n we locate bound states 
in the following way. For \m(k) > [Im(A;) < 0] f,nn{k,r) 
[f,n„(k,r)] converges as r —> °° because the right-hand side 
of Eq. (3) vanishes sufficiently quickly for the assumed long- 
range behavior of the potential [13]. The functions f^ n (k) = 
linv^oo '")lim(jt) < o are tne J° st functions familiar from 
general scattering theory [12, 13]. Introducing the matrices 
F ± (k) = \fnn(k)] we locate bound states by finding momenta 
k/, = i\kj, \ on the positive imaginary axis with vanishing deter- 
minant, i.e., 

det[F- {k b )\ = «■ ^f mn (h)c n (k h ) = 0. (5) 

n 

Iff the determinant of the Jost matrix F vanishes, then there 
exists a nonvanishing set of coefficients c„ fulfilling the right- 
hand side of Eq. (5). Moreover, at large distances the so- 
lutions u m (k[,,r) = Y*n§mn(kb,r)c„(kt,) vanish exponentially 
because u m (k b ,r) = 0[h+(k b r)} = 0[e-l**h > 0. There- 

fore, u m (kb,r) describes a bound state. Its binding energy is 
Eb = -\h\ 2 in units of h 2 /2/jA 2 . 

Now we are going to derive an approximate expres- 
sion for the binding energy in a weak anisotropic 2D 
potential V(r) = VqV(t), where Vq characterizes the 
strength of the potential. First, we introduce the ex- 
plicit expression for the regular solutions § m „(k,r) — 
jn{kr)b mn - Sq dr' g m {k,r,r')^ m 'V mm i{r)iSf m , n {k,r'), where 
g m (k,r,r') = k~ 1 [j m (kr)y m (kr J ) - j m (kr')y m (kr)} is the 
free Green's function of Eq. (1) [14]. Furthermore, the 
regular solutions of the Schrodinger Eq. (1) for k = are 
§mn(r) = r n+l l 2 h mn - fodr J g m (r,r J )Y im 'V mm ;{r)^ mln (r J ), 
where go{r,r') = \fr? ln(r' / r) and g m ^o(r,r') = 
(yrF /2m)[(r' / r) m — (r/r 1 )" 1 ] are corresponding free Green's 
functions. Now x /%(k/2)" +1 / 2 § mn (r)/n\ is the expansion of 
§ mn (k, r) to leading order in k. In Eq. (4) we replace § mn (k, r) 
and the scaled Hankel function by their respective leading 
order terms. Then we insert the asymptotic form of this ap- 
proximation for f mn {k) into the left-hand side of Eq. (5). The 
result is det[F" (Jfc)] = A [In (it) + y- ln(2) — in/2] + 1 — B = 0, 
where y is the Euler-Mascheroni constant. Solving this for k 



we find the binding energy 

E h = -4e- 2 ?- 2 



(6) 



We obtain the power series A = ^7=1-^ an d B — L/Li-B^ 
by expanding § mn (r) into a power series of the potential 
strength Vb, where j indicates the power of Vo. For a gen- 
eral (possibly anisotropic) 2D potential the leading terms of 
these series, including all partial waves, are given by 



A (1) = -/ drrVoo(r), 
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Here, i$(r) = S r dr>r»l 2 + m g n {r,r>)V nm {r>) + 

(1 /2n)r 1 / 2+ " fodr'r n+m -"V mn (r r ) and 
w m (r) = fodr , r J1 / 2 +">g m (r,r')V mn (r J ) 

Ln^o,m( 1 / 2n ) rl/2+ ' n Jo dr ' r ' v >m( r ')- We can calculate 
systematically higher order terms by expanding det[,F~ (k)] to 
higher orders in Vo. Very recently, Volosniev et al. [9] also 
found such an expansion of the binding energy by solving the 
2D Schrodinger equation directly for bound states. 

For concreteness let us now consider two polarized dipoles 
trapped in two parallel layers separated by a distance A. We 
define the z axis to be perpendicular and the x-y plane to be 
parallel to the layers. Without loss of generality, we assume 
that the dipoles are polarized within the x-z plane at an angle 
"& from the z axis. They interact via the interlayer potential 



V(r)=V - 



1 — 3(rcos9sini3' + cos'r}) 2 
(r 2 + 1)V2 



(11) 



Here, r = r(coscp,sincp) is the dimensionless projected vec- 
tor between the dipoles in polar coordinates (in units of A) 
and Vb = /Jd 2 /2%h 2 eoA is the interaction strength in units of 
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d/jt (rad) y (units of h 2 /2/jA 2 ) 



FIG. 1. (Color online) Binding energy of lowest interlayer bound 
states in the bilayer dipolar potential Eq. (11) as a function of (a) 
the polarization angle from the symmetry axis and (b) interaction 
strength Vo- The symbols mark numerical results, the dashed line 
an approximation in Ref. [1 1] for Vo = 1, and solid lines the corre- 
sponding approximation Eqs. (6)-(10). The latter approximation re- 
mains valid up to moderate Vo as long as \E/,\ <C 1. The insets show 
densities of the bound states at the indicated polarization angles and 
Vo = 10. The dark circles indicate the radius of one layer distance A. 

h 2 /2/jA 2 , with £o the electric constant and d the electric dipole 
moment (for magnetic dipoles Vo = /^W^ / '2%ti 2 A with /jq the 
magnetic constant and /jj the magnetic dipole moment). The 
potential fulfills AW = and / d 2 iV(r) = so that at least one 
bound state exists for all polarization angles. For $ = 0, this 
potential reduces to the central case discussed in Refs. [6, 7]. 
Using Eq. (6) and expanding A and B to fourth and second 
order, respectively, we recover the binding energy for perpen- 
dicular polarization given in Ref. [7]. 

In Fig. 1 we plot the binding energy of the lowest lying 
interlayer bound states in the potential (11) as a function of 
the polarization angle and interaction strength. Their bind- 
ing is strongest for perpendicular polarization and weakest 
for parallel polarization. We observe that the weak poten- 
tial approximation for the binding energy, Eqs. (6)-(10), re- 
mains valid up to moderate potential strengths Vo < 3 as long 
as the binding energy is sufficiently small, \Ejj\ <C 1. In con- 
trast, the approximation in Ref. [11], \Ej,\ ~ exp(l/c) with 
c = (l/87t 2 )/t/ 2 rt/Vv(r)ln|r-r'|V(r'), describes the angle 
dependence of the binding energy only poorly. Our numer- 
ical computations are based on Netlib's ZVODE solver, and 
we provide an explicit Jacobian for improved stability. We 
include partial waves up to sixth order. 

As an example, we consider bosonic 6 Li 40 K molecules sep- 
arated by A = 200nm. Then the energy scale in Fig. 1 is 
h 2 /2/jA 2 ~ 1.2 p-K (units of k B ) and V ~ 9.5 (cf. V = 10 in 
Fig. 1), with k B the Boltzmann constant. Therefore, the inter- 
layer bound state of LiK molecules should persist over a wide 
range of polarization angles under typical experimental tem- 
peratures in the nano-Kelvin regime. For fermionic 40 K 87 Rb 
molecules we have Vb ~ 1 with energy scale ~ 220 nK (units 
of kg). Therefore, it can be seen from Fig. 1(a) that they re- 
quire much lower temperatures. In contrast, for dipolar atoms 
with a magnetic dipole moment, Vb -C 1 (e.g., 52 Cr). Inter- 
layer dipoles with Vo -C 1 bind too weakly for all polarization 
angles to be stable at reasonable external parameters. 

The peak of the bound state wave function shifts along the x 



axis as the polarization direction changes from perpendicular 
to parallel [cf. insets in Fig. 1(a)]. This is because the attrac- 
tive dipole term in the DDI dominates over the quadrupole 
term. However, for $ > 0.4771, the quadrupole term domi- 
nates and two peaks develop. For parallel polarization, the 
dipole term vanishes and the two peaks become symmetric. 
In Fig. 1(a) the distance between these peaks is ~ 5A, typ- 
ically on the order of [im. If the stability requirements are 
fulfilled, this could simplify the creation of interlayer bound 
states at parallel polarization because of their greater overlap 
than states at perpendicular polarization. The parallel polar- 
ization bound state can be distinguished from the mainly per- 
pendicular one in a time of flight measurement. The time-of- 
flight image reflects the double peak in an asymmetric mo- 
mentum distribution. 



HI. SCATTERING 

Next we focus on the 2D scattering problem, for which 
k is real. Then the functions f^ n (k.r) attain the finite 
limit fm n (k) — lim,.^oo because the right-hand side 
of Eq. (3) vanishes at r — > °°. In order to capture the mix- 
ing of different partial waves, it is necessary to calculate the 
full S-matrix. We express the partial wave components of 
the scattering solution as a linear combination \)f mn (k,r) = 
L m ' §mm' (k, r ) c m'n (k) of the regular solutions, Eq. (2). Asymp- 
totically, we replace f , (k, r) in (j) mm / (k, r) with the Jost func- 
tions f~ m i(k). On the other hand, the general asymptotic scat- 
tering wave function in 2D is ^(k^) — > [e lk r + a(k,r) ^]/2tc, 

where fl(k,r) is the scattering amplitude. We match the two 
asymptotic expressions for the scattering wave functions by 
expanding the first exponential in ^(kjr) in terms of scaled 
Hankel functions. This way, we extract the S-matrix 

S(k)=F+(k)[F-(k)]- 1 (12) 

and the coefficients c mn {k) = [F~ {k)] m l- If the potential V(r) 
is central, S(k) is diagonal, with elements e 2l8m ^ and 5 m (k) 
the ra-th partial wave phase shift. 

Let us now derive approximate expressions for 
the phase shifts of anisotropic 2D scattering at low 
energies. First we introduce an iterative solution 
for the coefficients fm„(k,r) = LJ =0 fmn\k,r) as 
a power series in the potential strength Vo. From 
Eq. (4) and the expression for § ml ,(k,r) we obtain 

f^ J+l \k,r) = ±kfidr'h*(k r ')l*V^(T')$ n (k,r') 

and fim°\k,r) = 8 m „, with §mn(k,r) = 

- S^dSg^k, r, /) Zn< V W {'Wrtnik, r 1 ) and ^(k,r) = 
j„{kr)b mn . In the remainder of this section, we consider low 
scattering energies such that only up to two partial waves I 
and £' dominate the properties of the S-matrix. The phase 
shifts are given by tan28£(£) = ImSee/ReSce, where See is a 
diagonal matrix element of S and Sf is obtained by replacing 
£ — > £', By inserting Eq. (12) we obtain the phase shifts from 
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FIG. 2. (Color online) S-wave phase shift for anisotropic interlayer 
scattering of dipoles as a function of (a) the collision momentum 
for Vo = 1 (e.g., KRb molecules at A = 200 nm) and (b) the polar- 
ization angle. The symbols mark numerical results, solid lines the 
second-order Born approximation Eq. (14), dashed lines Eq. (15), 

and dotted lines Eq. (13) to second order in fmn ■ The weak po- 
tential approximation (13) describes scattering well from moderately 
small collision momenta and up to moderate potential strengths at all 
polarizations. The bound-state approximation (15) is only accurate 
at very low momenta and predominantly perpendicular polarization. 



the Jost matrix F since F + = F for real k. Then 



tan25 f (T)~2 



Rest 



(13) 



-U) 



with S(t = Su\6st(F )| 2 . Using terms f m ,Y' up to 
second order and expanding su to second order in 



Vo we find Im(s^) 



2/ (i) 7 (i) 
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jmjm ymjm' jut yut 1 ygJi 
= (l/k)Jodrp m (kr)V mn (r)q n (kr), 
( 1 A) Jo" dr P" ( kr ) tm V nm {r) Jq dr'g m (k, r, r')V mn (r J )j„ {kr') 
and p m , q n stand in for scaled Bessel functions of order m 
and n, respectively. The anisotropy of the potential enters the 

phase shift through mixing terms, such as Ij]) and the sum 



(21 

in Ip n . These terms are absent for a central potential so we 
recover a result in Ref. [6]. Furthermore, for small Vo we 
expand Eq. (13) to second order in Vo: 

tan28,(*) * -2/£> (*)+2[/f -®$ t -®^ (U) 

Thus, we recover a second-order Born approximation of the 
phase shift from our general formalism. We find a fur- 
ther approximation by considering the sum of all partial 
waves at small energies. For any S-matrix, cotJ^8^(fc) = 
Re(det5)/Im(det5) = -Re(detF~)/Im(detF~). As in the 
preceding section, we expand det,F~ around k = and use 
the fact that k is real to obtain 



cotVb((k) ~ - In 



k 2 
\E b \ 



(15) 



Here, Ej, is the binding energy, Eq. (6). 

In Fig. 2(a) we plot the s-wave phase shift of two dipoles 
with Vo = 1 interacting across two layers as a function of the 
collision momentum k, e.g., 40 K 87 Rb with different spins at 
A = 200 nm. The sharp increase of the s-wave phase shift at 
k — > (t3- = 0, 7t/6) is a consequence of the very weakly bound 
states. If So dominates in Eq. (15), a very small binding energy 
leads to a phase jump in So close to k — 0. Since the binding 
energy decreases strongly with increasing polarization angle, 
expression (15) describes the scattering of mainly in-plane po- 
larization only at unrealistically small collision energies. On 
the other hand, the weak potential approximation (13) de- 
scribes the numerics excellently at all considered momenta. 
It fails at large potential strengths Vo 3> 1 and small momenta. 
For large momenta k ^> 1 this approximation becomes identi- 
cal with the Born approximation (14). In Fig. 2(b) we observe 
that the s-wave phase shift can vary strongly with the polariza- 
tion angle at small momenta. The phase jump is caused by a 
weakly bound state at t3- ~ 0.37T, [see Fig. 1(a)]. This variability 
should be observable in the scattering of polar molecules. The 
bound-state approximation (15) describes this behavior qual- 
itatively even for such moderately large potential strengths as 
long as k <C 1 . The difference is mainly due to neglecting 
higher order partial waves in Eq. (15). For a small potential 
strength Vo <C 1 our approximation Eq. (13) describes scat- 
tering excellently at all polarizations. We find that this ap- 
proximation remains valid even at moderately large potential 
strengths Vo < 3 at higher energies k 2 > 1 . 



IV. CONCLUSIONS 

We have proposed a framework for calculating scattering 
and bound state properties for anisotropic 2D potentials. Our 
method generalizes the Jost formalism known from 3D scat- 
tering. We have derived systematic approximations for the 
scattering phase shifts and binding energy at low to moderate 
potential strengths. For weak potentials we have recovered a 
second-order Born approximation. The central equation (3) is 
also well-suited for numerical computations. 

We have applied our method to polar molecules trapped in 
a bilayer and polarized along an arbitrary direction. We find 
that absolute energies of 6 Li 40 K interlayer bound states are 
larger than their thermal energy in ultracold experiments even 
for nonperpendicular polarization. The s-wave phase shift of 
molecules with moderate or large DDI exhibits a strong de- 
pendence on the polarization angle. These results are im- 
portant, e.g., for the BEC-BCS crossover in fermionic polar 
molecules in bilayers [15, 16]. Varying the direction of the 
external polarizing field should influence the crossover from 
interlayer pair condensation to BCS pairing. 
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